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We theoretically study the Kelvin-Helmholtz instability in phase-separated two- component Bose- 
Einstein condensates using the Gross-Pitaevskii and Bogoliubov-de Gennes models. A flat interface 
between the two condensates is shown to deform into sawtooth or Stokes- type waves, leading to 
the formation of singly quantized vortices on the peaks and troughs of the waves. This scenario of 
interface instability in quantum fluids is quite different from that in classical fluids. 



I. INTRODUCTION 

A vortex sheet exists on the interface layer between 
two fluids with a relative velocity. When the relative ve- 
locity is sufficiently large, the vortex sheet becomes dy- 
namically unstable and the interface modes are amplified, 
which typically leads to roll-up patterns in the nonlinear 
stagei. This shear-flow instability is called the Kelvin- 
Helmholtz instability (KHI) and is related to several fa- 
miliar phenomena such as wind-generated ocean waves, 
flapping flags, billow clouds, and sand dunes. KHI is a 
fundamental instability in fluid dynamics and appears at 
various scales ranging from the laboratory to astronomi- 
cal scales, e.g. in plasmatic flows around the earth^ and 
relativistic flows in astrophysical jets^. 

Recently, KHI in quantum fluids, which we refer to as 
quantum KHI, has attracted growing interest. Quantum 
KHI was first realized by the Helsinki group^, in which 
quantized vortices penetrate from the A phase to the B 
phase of superfluid ^ He in a rotating cryostat. Quantum 
KHI is also a candidate mechanism of the pulsar glitches 
of rotating neutron stars with nucleon superflows^, the 
instability of the crystal-superfluid interface in ^He^, and 
the vortex formation in atomic Bose-Einstein conden- 
sates (BECs)i. A significant difference of quantum fluids 
from classical fluids is vortex quantization, and hence 
quantum KHI is expected to yield novel nonlinear dy- 
namics, which have not yet been explored. 

The quantum KHI observed in the Helsinki 
experiment^ can be explained by the hydrodynamic 
theory^. There exist two mechanisms that can cause 
quantum KHI. The first is dynamic instability (DI) 
induced by interface modes (ripplons) with complex 
frequencies, which is a direct analogue of the classical 
KHI. The second is the Landau (thermodynamic) 
instability (LI) caused by excitation of ripplons with 
negative energies (frequencies), which is absent in classi- 
cal KHI and thus unique to quantum KHI. We call these 
two mechanisms dynamic and thermodynamic KHI, 
respectively. In the helium experiment with dissipation 
caused by interaction with the normal component and 
the container wall, the thermodynamic KHI precedes 
the dynamic KHI because the critical velocity for the LI 



is generally lower than that for the DI^. In fact, the in- 
stability observed in Refj^ was actually thermodynamic 
KHI, in which the number of vortices was detected 
by NMR. However, the explicit nonlinear dynamics of 
vortex nucleation at the interface is unknown because of 
difficulties in directly observing vortex dynamics and in 
its hydrodynamic analysis. 

In this paper, we address the nonlinear dynamics by 
considering the shear-flow states in phase-separated two- 
component BECs of ultracold atoms. Atomic BECs are 
ideal systems to examine quantum KHI, since we can 
prepare an interface and shear-flow in a well-controlled 
manner, directly observe the vortex dynamics, and re- 
alize both dynamic KHI and thermodynamic KHI be- 
cause of the low dissipation rate. We employ the Gross- 
Pitaevskii (GP) and Bogoliubov-de Gennes (BdG) mod- 
els, which are more fundamental approaches than the 
hydrodynamic model in superfluid helium, for studying 
the stability of the shear-flow states and the vortex dy- 
namics. We flnd that the nonlinear dynamics of quantum 
KHI are quite different from those for classical KHI. Fig- 
ure [1] shows the typical dynamics of the quantum KHI. 
The initial state exhibits a flat vortex sheet on the in- 
terface between the two condensates [Fig. [I](a)]. When 
the relative velocity exceeds a critical value, the inter- 
face modes are amplifled to form sawtooth waves, with 
the vorticity being extremely localized at the peaks and 
troughs of the waves [Fig. [1] (c)]. Then, singly quan- 
tized vortices are released from the vortex sheet into the 
bulk of each condensate [Figs. [1] (d)-(f)]. Depending on 
the relative velocity and dissipation, we flnd a variety of 
nonlinear dynamics of quantum KHI. We also flnd that 
quantum KHI is observable in trapped BECs. 

II. SHEAR-FLOW STATE IN 
TWO-COMPONENT BOSE-EINSTEIN 
CONDENSATES 

We consider two-component BECs in a quasi- 
two-dimensional system under an external potential 
^j{y) U — l'^) for component j. In the mean-fleld 
theory at low temperatures, the two BECs are de- 
scribed by the condensate wave functions ^j(t,r) = 




FIG. 1: (color online) Nonlinear dynamics of the dynamic KHI in phase- separated two-component BECs without dissipation 
for the relative velocity Vr = V1—V2 = 0.98c > Vd. (Upper panels) The height and color show the vorticity cjeff and the density 
difference ni — 712 between the two condensates, respectively. (Lower panels) Two-dimensional plots of ni —712- (Right panels) 
The phases Gj (j = 1, 2) of the condensates for (f). The circular arrows show the rotational direction of the superflow around 
each vortex. 



^Jnj(t^ r)g*0i(^'r) ^\^]^ particle densities Uj and the 
phases 0j, which obey the coupled GP equations^ 



2m 



-v' + c/i + E^ifel^fcN (1) 



Here, rrij is the atomic mass and the coupling constants 
9iii 922^ and gi2 = 921 are related to the s-wave scat- 
tering lengths an, a22, and ai2 = a2i, respectively, by 
Qij — 2iT?i aij{mi + ruj) /ruimj {i^j = 1,2). In the sta- 
tionary state, we suppose that the two components un- 
dergo phase separation and form an interface layer along 
7/ = 0, sustained by the linear potential Uj = fji/ with a 
slight inclination fj . The condition for the phase separa- 
tion is given by /3 = 912/^/911922 > 1, and the densities 
ni and 77,2 are almost zero in the regions ^ > and ^ < 0, 
respectively, far from the layer. We consider here a suffi- 
ciently large value /3 = 10; the interface layer is thick and 
becomes indistinct when /3 ~ 1. Each component flows 
along the x-axis with the velocity Vj; the stationary wave 

function has a form = . / —ipj{y)e'^^^ , where ijjj is a 



real function and 6^ = 



-^VjX with the chem- 



ical potential fij. For simplicity, we set the parameters 
as /i = /ii - \miV^ = /i2 - 9 = 911 = 922, 

m = mi = m2, and / = /i = — /2- The units of length, 

\J I {vfi\i)^ T — Ti/iJi^ and 



time, and velocity are ^ 



c = ^y /i/m, respectively. 

Figure [2] (a) shows the profile of '^j(^) for / = 0.02/i/^. 
The shear-flow states are characterized by the effective 
superflow velocity Veff = (miniVi + m2n2V2)/(mini + 
m2n2) with Vj = ^V6j. The velocity field Veff(r) 
changes sharply at the interface from Vix {y < 0) to 
V2'k {y > 0). The effective vorticity cjeff = (V x Veff)^ 
is localized on the interface, which constitutes a vortex 
sheet. 
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FIG. 2: (color online) (a) Profiles of the wave functions of 
the two condensates (V^j) and the excitation (^V^j = Uj — Vj) 
in a shear-flow state. The real functions 6ipj for the lowest 
excitation with k — 0.29/(f for Vv — 0.79c are localized around 
the interface layer, (b) Dispersion relation ujq — uj — Ve^k of 
the lowest modes. The broken and solid curves show Im[c(;o {k)] 
for Vr/c = 0.98 and Re[a;o(A:)] for Vr/c = 0, 0.59, 0.79, 0.84, 
and 0.98, respectively, (c) Phase diagram of the quantum 
KHI for Vi — and V2 < 0. The curves are the boundaries 
of the DI and LI regions obtained by the BdG model (solid 
curves) and the hydrodynamic model (broken curves). 



III. LINEAR STABILITY ANALYSIS 
A. Hydrodynamic model 

First, it is instructive to discuss the stability of 
the shear-flow state with a hydrodynamic model. In 
this model, the interface layer is represented by a 
curve y = r]{t^x)^ neglecting its thickness. The La- 
grangian L = J dx J dy {Vi + 7^2 — ^12^1^2) of the two- 
component BECs is approximated with the interface- 
tension coefficient a and the interface area S as 



aS{r]), (2) 



3 



'A 

2 "^J 
2 



Ej with 



where Vj = f'^pt'^j - 

2^|V^jP + Uj{y)nj + ^djj'^'j' The interface tension 
is obtained by calculating the excess energy due to the 
presence of the interface^^. For small r^, we have S{r]) ~ 

J dx 1 + ^ {dxTff . Then the variation in Eq. (|2]) with 

respect to r] gives the equation Vi{r]) —V2{r]) -\- ad'^r] = 0, 
which corresponds to the Bernoulli equation in hydrody- 
namics. The motion of the interface can be determined 
by solving the Bernoulli equation with the kinematic 
boundary condition^. 

To understand the stability qualitatively, the profiles 
Gj — cx e~^~^^^^y cos{kx—ujt) of the phase fluctuations 
for the surface modes in a single-component BEC^^ are 
applied to each surface of the two condensates. As a 
result, for the interface modes with the form r] oc sm{kx — 
ujt), we obtain the dispersion relation. 
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where k > 0^ pj = mju'-^ and F = n[fi — 712/2 with 
the density near the interface in the Thomas-Fermi 
approximation rij ^ [p — Uj{0)]/g. The frequency u 
depends on the relative velocity Vr = Vi — V2 and the 
effective superflow velocity Vef[ = {piVi +^2^2)/ {pi + P2) 
at the interface of the stationary state. For 
our parameters we have pi = P2 = ^p/g^ 
F = Kff = (^1 + ^2)/2, and a = 

[2V2/3 - 0.514/3-1/^ - 2 (0.055/3-3/4 ^ 0.067/3-^/^) 
0.886^ according to RefJ^. 

For V^^ > 2^^^^^/Fa^ the imaginary part Im(cj) be- 
comes nonzero and the shear-flow states are dynami- 
cally unstable against excitation of the interface modes 
with k- < k < k-^ as in classical KHI, where k± = 

-V'^ . It is interesting 
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F lot with /co = 7w — T — Y 
to note that the dispersion relation of uj with V^^ 



has an inflection point as in the roton curve in superfluid 
helium, and Landau's argument of thermodynamic in- 
stability also applies to the present system. The effective 
velocity corresponding to the Landau critical velocity is 

-y? N P^ + p2- If V'eff exceeds 



given by 
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this value, there appear interface modes (ripplons) with 
negative energy quanta fiuj < {) and LI occurs in the pres- 
ence of dissipation. This instability is unique to quantum 
KHI since the shear-flow state cannot exist in thermal 
equilibrium in classical hydrodynamics. 

The interface modes with /c ^ 1/^ are not well de- 
scribed with the hydrodynamic model since the interface 
thickness ~ ^ is neglected there. In addition, for /c ^ 0, 
the phase velocity uj/k m Eq. (3) diverges to infinity. 
This unphysical divergence is caused by neglecting den- 
sity fluctuation of the condensates in this model, which 
is properly treated below. 



B. Bogoliubov-de Gennes model 

We can investigate the stability of the shear-flow 
states with the BdG model more precisely than the 
hydrodynamic model. The BdG model describes ex- 
citation modes in the bulk in addition to the inter- 
face modes, considering all degrees of freedom of the 
fluctuation. The fluctuation is written as a collec- 



tive excitation 5^ n 



around the shear- 



flow state Because of the translational symmetry 

along the x-axis, the excitation has a form (5^j(t,r) = 

giGO(t,x) \^^.(^y^^^-iut+ikx _^^^(^yyojt-ikx^ ^-^^Yi the norm 

V = ^jj(\uj\^ — \vj\^)dy > 0. The frequency 
cj(/c,Veff) = ^o(^) + Vei^k is Calculated by solving the 
BdG equations 
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Figure [2] (b) shows the numerical results for the dis- 
persion relations uJo{k) of the excitations with the lowest 
energies, which correspond to the interface modes since 
Uj and Vj are localized around the interface [Fig. [2] (a)]. 
In the present discussion, we restrict ourselves to the low- 
est mode, since the other modes have much larger criti- 
Hcal- -^efecities. A minimum in the curve Re(ct;o) appears 
for a finite and the curve reaches zero at K = ^d- 
When Vr > Vb^ the imaginary part Im(a;o) emerges for 
/cl < < /cr, where /cl (/cr) corresponds to k- (A:+) 
in the hydrodynamic model. This behavior is similar to 
that of the interface modes in the hydrodynamic model. 

For small k the lowest mode spreads over the whole 
system, and thus is not regarded as the interface mode. 
This mode causes the density fluctuation of the con- 
densates in the bulk, since the profile S'^j is continu- 
ously transformed to that of i/jj for A: ^ 0. The k = 
mode with zero energy always exists in the BdG modeP. 
Then, for /c ^ 0, the phase velocity ujo/k should be 
asymptotic to the phonon velocity, which is finite in the 
BdG model [see Fig. [2] (b)]. Since the phase velocity 
{uj-kVef^)/k (X 1/Vk in Eq. 1^ diverges, the BdG 
model is more accurate than the hydrodynamic model 
for k ^ 0. 

Figure [2] (c) shows the phase diagram for the DI with 
lm[uj{k)] ^ and the LI with uj{k) < for Vi = and 
V2 < 0, and then ]4ff = —Vr/2. The phase boundaries 
obtained from the BdG model (solid curves) agree with 
those from the hydrodynamic model with Eq. (|3j) (broken 
curves) except for k ^ 1/^ and k ^ 0. As increases, 
the LI region appears at K = and the boundary 
between the LI and stable regions is given by uj{k) = 0. 
The DI region is located above the LI region in Fig.[2](c). 
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The left and right boundaries between the DI and LI 
regions correspond to and kYi{Vr) in Fig. [2] (b) 

and meet at Vr = Vd- Note that the DI and LI regions 
never overlap. This is because the modes in the DI region 
which have zero norm u = and thus zero energy, do not 
cause energy dissipation and thus do not cause LI, where 
the energy of the fluctuations is formally written as huju. 



IV. NONLINEAR DYNAMICS 

When Vr > Vd or Vr > Vl^ the interface modes are 
amplified, leading to nonlinear dynamics. Note that the 
thermodynamic KHI triggered by LI is not described by 
the energy-conserving GP Eq. ([T]). Here, we introduce a 
dissipative model, which can describe the nonlinear dy- 
namics triggered by both DI and LI^^. The dissipative 
nonlinear dynamics is obtained by solving the modified 
GP equation 



2mj 



where 7 > is the phenomenological dissipation 
constant^^. In the linear stage of the instability, 
the DI and LI modes are exponentially amplified as 
exp[|Im(a;)|t] and exp[7|cj|t], respectively. Thus the dy- 
namics crucially depend on 7 and Vr. 



A. Dynamic KHI 

We now demonstrate the dynamic KHI with 7 = for 
Vr = 0.98c > Vd, where LI is absent because of the ab- 
sence of dissipation. We numerically solve the GP equa- 
tion in a periodic system with a period L = 64^ along the 
X-axis. The initial state is with a small random seed 
to trigger the instability [Fig. [1] (a)]. In the linear stage 
of the instability, the sine wave with ~ 0.29, which has 
the maximum imaginary part max/g {Im [(^(/c)]}, is dom- 
inantly amplified. As the amplitude becomes large, the 
sine wave is distorted by nonlinearity [Fig. [1] (b)], and 
deforms into a sawtooth wave [Fig. [1] (c)]. We note that 
the vorticity cjeff increases on the edges of the sawtooth 
waves and creates singular peaks [Fig. [1] (d)]. Subse- 
quently, each singular peak is released into each bulk, 
becoming a singly quantized vortex with a circulation of 
hi = h/m [Fig. [11(e)]. 

The release of vortices reduces the vorticity of the 
vortex sheet and therefore reduces the relative velocity 
across the interface. The relative velocity can be roughly 
estimated with the total vorticity of the vortex sheet di- 
vided by the length of the sheet. Since six singly quan- 
tized vortices are released from the edges of the sawtooth 
waves [Fig.[T](f)], the relative velocity decreases by about 
6k,/ L ^ 0.6c to below the threshold Vb for DI, so that 
additional vortex creation is suppressed. The released 
vortices drift along the interface and the system never 



recovers the initial fiat interface. These nonlinear dy- 
namics are quite different from those in classical KHI, 
where the interface wave grows into roll- up patterns^-. 



B. Thermodynamic KHI 

We next consider the thermodynamic KHI with dissi- 
pation 7 = 0.03 for Vl < Vr = 0.79c < Vb [see Fig. [3] 
(top)], where DI is absent. In the nonlinear stage, we 
find that the interface has fiattened troughs and peaked 
crests [Fig. [3] (a)] like the Stokes wave known as a finite- 
amplitude wave at the surface of water^. In this case, 
the vorticity is highly localized at the crests [Fig. [3] (b)] 
and a singly-quantized vortex is released from each crest 
into the upper side [Fig. [3](c)], in contrast to the case of 
the dynamic KHI where the vortices are nucleated in the 
both sides. 

The dissipation drags the released vortices away from 
the interface to slow down the upper phase [Fig. [3](d)]. 
Such a slowdown due to the release of quantized vortices 
is interpreted as phase slippag As a result, four singly 
quantized vortices are released, which reduces the relative 
velocity across the interface by about Ak/L ~ 0.4c to 
below the threshold Vi, for LI. Then the vortex nucleation 
stops and the vortex sheet recovers a uniform line with 
less vorticity than that of the initial vortex sheet. 



(a)//r=844 {b)t/T=9U (c) t/T=93S (d) ^/2-1250 




FIG. 3: (color online) Dynamics of the density difference 
ni — 712 in the quantum KHI with dissipation 7 = 0.03 for 
Vi < K = 0.79c < (top) and Vr = 0.98c > (bottom). 
The field of view is 32^ x 64^. 



C. Coexistence of dynamic and thermodynamic 
KHIs 

For Vr = 0.98c > T^d with dissipation 7 = 0.03, the 
dynamic and thermodynamic KHIs coexist. Figure [3] 
(bottom) shows typical dynamics after the linear devel- 
opment, in which the growth rate |Im(co')| of the DI is 
much larger than the rate 7|cj| of the LI. In the early 
stage [Fig. [3] (e)], the dominant DI causes modulation of 
the interface as in Fig.[T](b), but the dissipation deforms 
the sawtooth wave [Fig. [3] (f)]. 

Depending on the wave number of the dominantly am- 
plified mode, three vortices are released upward from the 
peaks of the wave as shown in Fig.[3](g). As in Figs. [3] (a)- 
(d), the vortices released from the interface move upward 
due to the dissipation. In this case, the relative velocity 
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across the interface decreases by ~ 0.3c between 

the thresholds and Vl. Then the dynamic instabihty 
stops and the thermodynamic instabihty becomes dom- 
inant, forming the Stokes-type waves [Fig. [3] (h)]. The 
subsequent dynamics are similar to that in Figs. [3] (b)- 

(d). 



D. Quantum KHI in trapped condensates 



where the component 2 with 30 vortices surrounds the 
component 1 with no vortices in the initial state. Such 
an initial state may be prepared by using a Laguerre- 
Gaussian beam with orbital angular momentum^^ or the 
vortex pump proposed in Ref.^^. The role of the force F 
in Eq. (j3j) is played by the centrifugal force for the rotat- 
ing component 2. We see that the dynamic KHI appears 
on the circular interface. 




FIG. 4: (color online) Dynamics of the density difference 
ni — 712 in the quantum KHI without dissipation in a trapped 
system. The total number of atoms is 3.1 x 10^ with an equal 
population in each component. The trap potential is Uj{Y) = 
f (cj^r^ + ujzz^) with = + 2/^, cjx = 27r x 80Hz, and 
ujz — 2n X 4kHz, and the simulation is performed in quasi- two 
dimensions. The field of view is 72/im x 72/im. 

Experimentally, quantum KHI can be realized in a 
trapped system. We assume that the components 1 and 
2 are the hyperfine states \F = I^tuf = —1) and |2, 1) of 
^''Rb atoms, where an : ai2 : a22 = 0.97 : 1 : 1.03 with 
their average being 5.5 nm according to Refs.^^^^^. Fig- 
ure [4] shows the dynamics of the condensates obtained by 
numerically solving the GP equation without dissipation. 



V. CONCLUSIONS 

We have studied quantum KHI in phase-separated two- 
component BECs with shear-flows using the GP and BdG 
models and found that a variety of interface patterns are 
formed, such as sawtooth and Stokes-type waves, and 
that quantized vortices make the nonlinear dynamics dis- 
tinct from those in classical fluids. The study of quan- 
tum KHI in three-dimensional systems is an interesting 
direction for future work, where the vortex reconnections 
would change the nonlinear dynamics qualitatively from 
those in two-dimensional systems. Details of quantum 
KHI in trapped BECs will be published elsewhere. We 
expect that various instabilities in classical fluid mechan- 
ics may have counterparts in quantum fluids, in which 
quantized vortices could yield novel nonlinear dynamics. 
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